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Abstract: The neutrino quasinormal modes of the Reissner-Nordstrom (RN) black hole are 
investigated using continued fraction approach. We find, for large angular quantum number, that 
the quasinormal frequencies become evenly spaced and the spacing of the real part depends on the 
charge of the black hole and that of the imaginary part is zero. We then find that the quasinormal 
frequencies in the complex lo plane move counterclockwise as the charge increases. They get a 
spiral-like shape, moving out of their Schwarzschild value and “looping in” towards some limiting 
frequency as the charge tends to the extremal value. The number of the spirals increases as the 
overtone number increases but it decreases as the angular quantum number increases. We also find 
that both the real and imaginary parts are oscillatory functions of the charge, and the oscillation 
becomes faster as the overtone number increases but it becomes slower as the angular quantum 
number increases. 
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1. introduction 


It is well-known that the quasinormal frequencies of a black hole are defined as proper solutions of 
the perturbation equations belonging to certain complex characteristic frequencies which satisfy 
the boundary conditions appropriate for purely ingoing waves at the event horizon and purely 
outgoing waves at infinity jlj] . The quasinormal modes (QNMs) are entirely fixed by the structure 
of the background spacetime and are irrelevant of the initial perturbations [1]. Thus, it is gen¬ 
erally believed that QNMs carry a unique footprint to directly identify the existence of a black 
hole. Through the QNMs, one can extract information about the physical parameters of the black 
hole—mass, electric charge, and angular momentum—from the gravitational wave signal by fitting 
the few lowest observed quasinormal frequencies to those calculated from the perturbation anal¬ 
ysis. Meanwhile, the study of QNMs may lead to a deeper understanding of the thermodynamic 
properties of black holes in loop quantum gravity [||] [)|, as well as the QNMs of anti-de Sitter 
black holes have a direction interpretation in terms of the dual conformal field theory i § ®- 
Therefore, much attention has been devoted to the study of the quasinormal modes in the recent 
thirty years f7f-f2l||. 

Although QNMs of scalar, electromagnetic and gravitational perturbations are studied exten¬ 
sively, the investigation of the QNMs of the Dirac field is very limited |22| p3| [ 24 1. Furthermore, 
in these papers the study of the Dirac QNMs were limited to use WKB or Poshl-Teller potential 
approximative methods. The WKB and Poshl-Teller potential approximative methods may be 
valuable tools in certain situations, but in general results obtained from the approximative meth¬ 
ods should not be expected to be accurate. It is well known that the Leaver’s continued fraction 
technique is the best workhorse method to compute highly damped modes and is very reliable. 
However, the method can not be directly used to study the QNMs of the Dirac fields if we take the 
standard wave equations. The reason to obstruct the study of the Dirac QNMs using continued 
fraction technique is that for static black holes the standard wave equation 


( d 2 


Z± = K±Z± 


( 1 . 1 ) 
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possesses a special potential 


2 a 

v± = A 2 —r ± A 


d y/A 
dr* r 2 


( 1 . 2 ) 


which is the function of y/A ^2|[^4], where A is a function related to the metric, say A = 
r 2 — 2 Mr + Q 2 for the RN black hole. We can not investigate the Dirac QNMs using the continued 
fraction approaches because we have to expand the potential as a series at the event horizon but 
the factor y/A does not permit us to do that. 


Recently, we [25] found that the wave function and potential of the Dirac field can be expressed 
as new forms, and the new wave function are related to the standard one in a simple way [p5|[. 


Starting from the new wave function and potential, we 25, 26] studied that the Dirac QNMs of the 
Schwarzschild black hole using continued fraction [27] and Hill-determinant approaches |2S|ll2^ 


and 

the Dirac QNMs of the Schwarzschild-anti-de Sitter and Reissner-Nordstrom-anti-de Sitter black 
holes using Horowitz-Hubeny approach |H|. Considering the RN black hole spacetime provides 
a better background than the Schwarzschild geometry, in this paper we will extend the study in 
Ref. (2f| to the RN black hole, and to see how the quasinormal frequencies depend on the charge 
of the black hole. 

The organization of this paper is as follows. In Sec. 2 the decoupled neutrino equations and 
corresponding wave equation in the RN spacetime are obtained using Newman-Penrose formalism. 
In Sec. 3 the numerical approach to compute the neutrino QNMs is introduced. In Sec. 4 the 
numerical results for the neutrino QNMs in the RN black hole are presented. The last section is 
devoted to a summary. 


2. Neutrino equations in the Reissner-Nordstrom spacetime 


The neutrino equations [3C] can be expressed as 


y/2 V BB 'P b = 0, 

V2Vbb'Q B = 0 , 


( 2 . 1 ) 


where V bb 1 is covariant differentiation, P B and Q B are the two-component spinors representing 
the wave functions, Pb' is the complex conjugate of Pb■ In the Newman-Penrose formalism |3l]] 
the equations become 


(D + e- p)P° + (6 + vr - a)P 1 = 0, 

(A + n - 7 ) P 1 + (5 + (3- t)P° = 0, 

( D + e — p)Q 0 + (6 + if — a)Q^ = 0, 
(A + jl — 7 )Q^ + {5 + f3 — f)Q ^ = 0. 


The null tetrad for the RN black hole can be taken as 

^ = 1 , 0 , 0 ), 

„ 1, A 

n^ = - 1 , 0 , 0 

2 r z 


1 


y/2r 


0 , 0 , 1 ,- 


sin6 


( 2 . 2 ) 


(2.3) 
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with 


A = r 2 — 2 Mr + Q 2 = (r — r + )(r — r_), (2-4) 

where M and Q are the mass and charge of the RN black hole, r + = M + y/ M 2 — Q 2 and 
r_ = M — \JM 2 — Q~ represent the event and inner horizons. The wave functions can be expressed 


as 


P° = -K_ 1/2 (r)S_ 1/2 (0)e 
P 1 =R +1/2 (r)S +1/2 (9)e 
Q 1 = ^+l/2( r )*S'-l/2(^) e 


Q°' = —M _1 /2 (r)S +1/2 (9)e 




—i{ujt—rrup) 


—i(ut—rrup) 


—i(ut—mp) 


(2.5) 


where co and m are the energy and angular momentum of the neutrino particle. After the tedious 
calculation Eq. (2.2) can be simplified as 


with 


yfK.T>Q^_i/2 = AVAM + 1 / 2 , 

( 2 . 6 ) 

^/ADq ( \/AR+ 1/2 ) = AM_x/ 2 ) 

(2.7) 

^l/ 2 'S’+l /2 = — ^5'-l/2) 

( 2 . 8 ) 

^i/ 2 ^- 1/2 = ^S+i/ 2 , 

d ir 2 co n dA 

n dr A A dr ’ 

(2.9) 

1 d ir 2 to n dA 
n dr A A dr ’ 

dm 

£n = ao + . „ +ncotfc», 

dd sm (J 

„+ d m 

£ n= ns, ■ a + ncot0 ’ 

dd sm (J 

( 2 . 10 ) 


where A = ± (l + (l is the quantum number characterizing the angular distribution). We will 
focus our attention on the massless Dirac field in this paper. Therefore, we can eliminate M_i / 2 
(or \/AR + 1 / 2 ) from Eqs. (^) and fl2.7|) to obtain a radial decoupled Dirac equation for \/AK+i /2 
(or M_ 1 / 2 )- Then, introducing an usual tortoise coordinate 

'2 


dr* = A^ r ’ 


and resolving the equation in the form 


A~ s / 2 


we obtain the wave equation 


d 2 * s 

dr 2 


-4' 


+ (w 2 — V r s )4' s — 0, 


with 


4r 2 dr 


,d_ 

dr 


+ 


s 2 r 4 


d_ 

dr 


2 d 

+ isuir — 
dr 


—7 + 


a 2 a 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 


We will study the neutrino quasinormal frequencies of the RN black hole using Eqs. (2.12) and 
(2.14) and with the help of the continued fraction approach. 
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3. Numerical Approaches 


For the RN black hole, the QNMs are defined to be modes with purely ingoing waves at the event 
horizon and purely outgoing waves at infinity jlj] . Then, the boundary conditions on wave function 
at the horizon (r = r + ) and infinity (r —> +oo) can be expressed as 




(r — r .|_) 2 2k 

y—s+iu gicor 


r -> r + , 
r —> Too, 


(3.1) 


where = 


r|—r 


2 ri 


is the surface gravity on the event horizon r + . 


A solution to Eq. (|2.12 ) that has the desired behavior at the boundary can be written in the 


form 


= r(r - r_|_) 


'f-, .-l-f+2 iu+ 

+ (r — r_) 


2 K_ e iu{r—r-) 


E' 

m=0 


r — r+ 


(3.2) 


where = 1 , J is the surface gravity on the inner horizon t~_. If we take r + +r_ = 1 and 

b = r. |_ — r_, the sequence of the expansion coefficients { a m : m = 1,2,....} is determined by a 
three-term recurrence relation staring with oq = 1 : 


aoai T Po&o — 0, 

T Pm®m T 7m®m-1 = 0) ?71 = 1, 2, — (3-3) 

The recurrence coefficient a m , p m and are given in terms of m and the black hole parameters 
by 


oi m = T (Co T 1)7 tz T Co, 

Pm = —2m 2 T (Ci T 2)771 T C 3 , 

7 m = T - (C 2 — 3)777 T C 4 — C 2 T 2, (3-4) 


and the intermediate constants C m are defined by 

iu(r\ T T’ 2 ) 


Cq = 1 — s — iuj — 


Ci = —4 T 2iuj(2 T b) + 


2 iu(r\ T r 2 ) 


C 2 — s T 3 — 3io; — 


7 w(ri T r 


C 3 = or (4 T 26 — 4 r + 7 ’_) — s — 1 T (2 T 6 )io; — A 2 T (2u; T i ) 


w(r^_ T T’ 2 ) 


C 4 = s T 1 — 2a ; 2 — (2s T 3)?'o; — (2a; T i ) 


a(r 2 Tr 2 ) 


(3.5) 


It is interesting to note that the three-term recursion relation is obtained form the wave equation 
(2.13) directly. The calculation is simpler than other cases. For example, the coefficients of the 
expansion for the electromagnetic and gravitational fields in the RN black hole are determined by 
a four-term recursion relation which should be reduced to a three-term relation using a Gaussian 
eliminated step (see [pi for details). 
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The series in (3.2) converges and the r = +oo boundary condition ®3) is satisfied if, for 
given s and A, the frequency lo is a root of the continued fraction equation 



1'1/m Otm—2'Ym—l C^OTl 

Pm —1 Pm—2 Po 


Q^mTm+l O t m+l n /m+2 ®m+2'7m+3 

Pm+1 Pm+2 Pm +3 


(m = 1,2...). (3.6) 


This leads to a simple method to find quasinormal frequencies— defining a function which returns 
the value of the continued fraction for an initial guess at the frequency, and then use a root finding 
routine to find the zeros of this function in the complex plane. The frequency for which happens 
is a quasinormal frequency. The nth quasinormal frequency is usually found to be the most stable 
root of the nth inversion |27t|. 


4. Neutrino Quasinormal Modes of the Reissner-Nordstrom black hole 

In this section we will present the numerical results of the neutrino quasinormal frequencies of 
the RN black hole obtained by using the numerical approach just outlined in the previous section. 
The results will be organized into three subsections: the fundamental quasinormal frequencies, the 
QNMs for large angular quantum number and highly damped modes. 

4.1 Fundamental quasinormal frequencies 

We find that the change of the neutrino quasinormal frequencies becomes faster as the charge 
increases. Therefore, as we approach the extremal value Q = 1/2, the convergence of the continued 
fraction method slower, and the required computing time gets longer. The fundamental neutrino 
quasinormal frequencies of the RN black hole for Q = 0 to Q = 0.499 and A = 1 are listed in 
the table I. We know from the table that the real part of the fundamental quasinormal frequency 
increases as Q increases, but the imaginary part decreases first and then increases as Q increases 
(we should note that the magnitude of the imaginary part increases first and then decreases). It 
is shown that the intermediate decay of the Dirac perturbation around the RN black hole depends 
on the charge Q. The oscillating frequency of this decay which corresponds to ojr increases as Q 
increases, and the amplitude per unit time ~ e^ 7 increases first and then decreases. We should 
point out that this property is true of all fundamental quasinormal frequencies with different 
angular quantum number. 

4.2 QNMs for large angular quantum number 

The quasinormal frequencies for A = 1 to A = 40 and n = 0, 1, 2 are described by Fig. [I] which 
shows that A uj = cua+i — co\ as a function of A. From the figure we know that, for the cases of 
n = 0, n = 1 and n = 2, the quasinormal frequencies become evenly spaced for large A and the 
spacing is given by 

A u(Q = 0.0) = 0.38490 - 0.0000b 
Aw(Q = 0.1) = 0.38750 - 0.0000b 
A co(Q = 0.2) = 0.39581 - 0.0000b 

A w(Q = 0.3) = 0.41170 - O.OOOOi. (4.1) 

That is to say, the spacing of the real part is related to the charge Q and the spacing of the 
imaginary part becomes zero for large A. 
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Q 

<jj 

Q 

Ul 

Q 

UJ 

0.000 

0.365926 -0.193965i 

0.410 

0.427193 -0.196008i 

0.460 

0.452381 -0.190363i 

0.050 

0.366591 -0.194066i 

0.415 

0.429341 -0.195734i 

0.465 

0.455387 -0.189224i 

0.100 

0.368615 -0.194366i 

0.420 

0.431562 -0.195418i 

0.470 

0.458469 -0.187904i 

0.150 

0.372096 -0.194851i 

0.425 

0.433862 -0.195052i 

0.475 

0.461608 -0.186375i 

0.200 

0.377208 -0.195488i 

0.430 

0.436242 -0.19463i 

0.480 

0.464769 -0.184605i 

0.250 

0.384241 -0.196213i 

0.435 

0.438707 -0.194144i 

0.485 

0.467892 -0.18257H 

0.300 

0.393651 -0.196891i 

0.440 

0.441258 -0.193585i 

0.490 

0.47089 -0.180287i 

0.350 

0.406184 -0.19722i 

0.445 

0.4439 -0.192942i 

0.495 

0.473689 -0.17785H 

0.400 

0.423108 -0.196443i 

0.450 

0.446635 -0.192202i 

0.497 

0.474767 -0.176866i 

0.405 

0.425116 -0.196243i 

0.455 

0.449462 -0.191348i 

0.499 

0.475833 -0.175873i 


Table 1: The fundamental neutrino quasinormal frequencies of the RN black hole for Q = 0 to Q = 0.499 
and A = 1. As charge Q increases, the real part of the fundamental quasinormal frequency increases, but 
the imaginary part decreases first and then increases. 



Figure 1 : The spacing Aw = oj\ + i — uj\ as the functions of A for n = 0, 1 , 2 quasinormal frequencies. 
The left three figures are drawn for A u>R e = Re(ui\+i) — Re(u>\ ) which show that the spacing of the real 
part is related to the charge Q for large A. The right three figures are for A= Jm(w a+i) — Im{u>\) 
which show that the spacing of the imaginary part becomes zero for large A. 


4.3 Highly damped modes 

The Fig. [2] is drawn for n = 0 to n = 3 and A = 1, and Fig. |j is for n = 5,6,7,8,10 and 
A = 2. In these figures the quasinormal frequency trajectories are formed only by points which are 
quasinormal frequencies at increments in Q of 0.01 for Q = 0 to Q ~ 0.40 and at increments in Q 
of 0.001 for Q ~ 0.40 to Q = 0.495. We now proceed to discussed the figures in more detail. 
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The left columns in the Figs. || and || describe the behavior of quasinormal frequencies in 
the complex ui plane which show that the frequencies generally move counterclockwise as the 
charge increases. They get a spiral-like shape, moving out of their Schwarzschild value (Q = 0) 
and “looping in” towards some limiting frequency as the charge Q tends to the extremal value, 
Q = 1/2. For a given A, we observe that the number of spirals increases as the overtone number 
increases. However, for a given overtone number n, increasing A has the effect of “unwinding” the 
spirals, as we see in the two figures that the spiral begins at n = 2 for A = 1 but it starts at n = 6 
for A = 2. 

The second and last columns in the Figs. |2| and |3| illustrate that the real and imaginary parts 
of the quasinormal frequencies are the functions of the charge Q. We know from these figures that 
both the real and imaginary parts of the frequencies are oscillatory functions of the charge. The 
oscillation starts earlier and earlier as the overtone number n grows for a fixed A, but it begins 
later and later as the angular quantum number increases for a fixed n. Meanwhile, the oscillation 
becomes faster as the overtone number increases for a given A, but it becomes slower as A increases 
for a given n. 


5. summary 

The wave equation for the neutrino fields in the RN black hole spacetime is obtained by means 
of the Newman-Penrose formulism. The expansion coefficients of the wave equation satisfying 
appropriate boundary conditions are determined by a three-term recurrence relation. Then, the 
neutrino quasinormal frequencies of the RN black hole spacetime are evaluated using continued 
fraction approach and the results are presented by table and figures. We find, for any angular 
quantum number, that the real part of the fundamental quasinormal frequency increases as the 
charge increases, but the imaginary part decreases first and then increases as the charge increases. 
The fact tell us that the intermediate decay of the Dirac perturbation around the RN black hole 
depends on the charge. As the charge increases, the oscillating frequency of this decay increases and 
the amplitude per unit time increases first and then decreases. We then show that the quasinormal 
frequencies become evenly spaced for large A and the spacing of the real part is related to the charge 
and that of the imaginary part becomes zero. We also find that the frequencies in the complex 
Lo plane generally move counterclockwise as the charge is increased. They get a spiral-like shape, 
moving out of their Schwarzschild value and “looping in” towards some limiting frequency as the 
charge tends to the extremal value. For a given A, we observe that the number of spirals increases 
as the overtone number increases. However, for a given overtone number, increasing A has the effect 
of “unwinding” the spirals. At last, we find that the real and imaginary parts of the quasinormal 
frequencies are oscillatory functions of charge, and the oscillating behavior starts earlier and earlier 
as the overtone number grows for a fixed A, but it begins later and later as A increases for a fixed 
overtone number. 
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Figure 2: The left 4 figures describe the behavior of quasinormal frequencies for n = 0 to n = 3 and A = 1 
in the complex w plane which show that the frequencies generally move counterclockwise as the charge is 
increased and the number of spirals increases as the overtone number increases. The others graph Re(io) 
and Im(iS) of the quasinormal frequencies versus charge Q which tell us that both the real and imaginary 
parts are oscillatory functions of charge and the oscillations become faster as the overtone number increases. 
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Figure 3: The left 5 figures describe the behavior of quasinormal frequencies for n = 5,6,7,8,10 and 
A = 2 in the complex to plane and the others graph Re(u ) and Im{ui) of the quasinormal frequencies versus 
charge Q. By comparing it with Fig. || we know that the number of the spirals decreases and the oscillations 
become slower as the angular quantum number increases. 
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